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Gluonic pole matrix elements explain the appearance of single spin asymmetries (SSA) in high- 
energy scattering processes. They involve a combination of operators which are odd under time 
reversal (T-odd). Such matrix elements appear in principle both for parton distribution functions 
and parton fragmentation functions. We show that for parton fragmentation functions these gluonic 
pole matrix elements vanish as a consequence of the analytic structure of scattering amplitudes in 
Quantum Chromodynamics. This result is important in the study of the universality of transverse 
momentum dependent (TMD) fragmentation functions. 
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Cross sections for high energy scattering processes are 
given by a convolution of partonic cross sections with 
parton distribution functions (PDF) and parton frag- 
mentation functions (PFF). These functions interpreted 
as momentum distributions and parton decay functions, 
respectively, are given as matrix elements of quark and 
gluon operators [l|-|4[. Among such observables, gluonic 
pole or Qiu-Sterman matrix elements involve a combi- 
nation of operators which is odd under time reversal 
(T-odd). These matrix elements have been extensively 
studied [oT-[l^|. They explain the appearance of single 
spin asymmetries (SSA) in high-energy scattering pro- 
cesses. These correlation functions show up in combi- 
nation with calculable hard parts that may differ from 
the partonic cross sections through specific calculable fac- 
tors and signs @, [H, 14|. In this paper we use general 
properties of scattering amplitudes in Quantum Chro- 
modynamics (QCD) to study the support properties of 
these parton correlation functions. Specifically, assuming 
unitarity and analyticity properties to hold for forward 
parton-hadron scattering amplitudes, we uncover their 
singularity structure which then determines the corre- 
sponding properties of the PDFs and PFFs. Using this 
analysis for fragmentation functions, we show that sin- 
gle gluon and multi-gluon pole matrix elements vanish 
in the limit when the momenta of these gluons become 
zero. Since these multi-gluonic pole matrix elements ap- 
pear in integrated and weighted transverse momentum 
dependent (TMD) fragmentation functions, as a conse- 
quence of them being zero, all leading T-odd effects in 
the matrix elements are all part of the final state interac- 
tions among the fragmentation remnant and final state 
hadron 0, [l5| rather than from T-odd partonic operator 
combinations. Thus, when transverse momentum depen- 
dent (TMD) T-odd fragmentation functions appear in 
observables like SSA they are convoluted with the stan- 
dard partonic cross sections. The vanishing of gluonic 



pole matrix elements for any number of gluons, provides 
a general proof of the universality of these TMD frag- 
mentation functions [Tfj| | . 

Spectral studies of the gluonic pole matrix elements for 
fragmentation, specifically for quark-qua rk-g luon matrix 
elements with just one gluon field [17|, uM, already in- 
dicated that they vanish. Our arguments presented here 
are more general and can be applied to multi-gluonic pole 
matrix elements. They depend only on the analytic struc- 
ture of scattering amplitudes, yet they do not depend on 
the details of partonic or hadronic masses, and are insen- 
sitive to integrations over transverse momentum. 

We begin our analysis by considering high-energy scat- 
tering processes where the structure of hadrons is ac- 
counted for using quark and gluon correlators, which are 
Fourier transforms of forward matrix elements of non- 
local quark and gluon operators between hadronic states. 
For instance the quark-quark correlator 

x(P|^(0)W M ^(0|P>J LP , (i) 
where = Vcxp[~ ig J c ds- A a (s) t a ] is the gauge link 



The non- locality 



that ensures gauge invariance 
of the matrix elements needed to describe the distribution 
functions is restricted to the light-front (LF) and it is 
convenient to use the Sudakov decomposition 



k = xP + an + k T , 



(2) 



in terms of a generic light-like four-vector n satisfying 
n 2 = and P ■ n = 1. In a particular hard pro- 
cess, its role is played by other momenta that are hard 
with respect to the hadron under consideration, e.g. 
?i ps P' I P-P'. We can then also work with light-cone co- 
ordinates. Including mass effects one would have n_ = n 



and n + = P 



\M 2 n; with fc ± 



k ■ 
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FIG. 1: The quark-quark correlators that establish the 
non-perturbative connection between partons and hadrons for 
distribution functions (a) or fragmentation functions (b) and 
a quark-quark-gluon (multi-parton) correlator (c). 



k + = k-n = x and k" = k ■ P - \ xM 2 = a + \ xM 2 . 
The transverse momentum is orthogonal to n and P. In a 
hard process the dependence on k~ of a particular corre- 
lator is not important and it is integrated over, leaving us 
with the restricted light-front non- locality £ + = (LF). 
The expansion of these correlators (in Dirac space) con- 
tains the TMD distribution functions depending on mo- 
mentum fraction x and transverse momentum k 2 . Upon 
integration over k T one obtains the collinear correlators 



2tt 



Axf-P 



p\mu$. 4] m\p)\ LC > (3) 



where non- locality is restricted to the light-cone (LC: £ ■ 
n = £ T = 0) and the gauge link is unique, being the 
straight-line path along n. These collinear correlators are 
expanded in the 'standard' parton distribution functions, 
depending solely on the momentum fraction x. Here we 
will not discuss the scale dependence 23-261]. 

The quark-quark light-front correlator that plays a role 
in the fragmentation of partons is 



^(z,k T ) 



E 

X 



d(t-P)d 2 £ T 
(2tt) 3 
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with the quark momentum, k = \ P + k T + a n, i.e. a 
Sudakov expansion with x = 1/z > 1. In this case one 
often refers to the hadron transverse momentum P± = 
~zk T (in a frame in which the parton does not have 
a transverse momentum (fcj^ = 0)). Diagrammatically 
these correlators are represented in Fig.[T] In this paper, 
we start the investigation of multi-parton correlators by 
looking at the case with one additional gluon, as given 
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FIG. 2: Integrating parton correlators over k~ allows connect- 
ing them to a single anti-parton - hadron scattering four-point 
function A(k 2 ;s,u) (middle). Depending on the value of x, 
the imaginary part of this amplitude represents the (anti)- 
parton distribution or fragmentation correlators. 

in Fig. [1] (c). They appear in azimuthal asymmetries 
involving the fc r -weighted correlator, 

1>q [U] (x) = [ d 2 k T k«& u \x,k T ) 



(5) 



= ig(x) + C^7r*g(x,x) 



which is decomposed in pieces <&g and $g, that con- 
tain T-even and T-odd operator combinations, respec- 
tively [22j. The T-odd parts come with calculable glu- 
onic pole factors Cq^ that depend on the gauge link. In 
this paper we focus on the connection of this part to the 
zero momentum [x\ — > 0) limit of a quark-quark-gluon 
correlator 



dtf-P) d(r,-P) 



2tt 2tt 

x (P\mu^ v] gG^^U^m |P)J LC , (6) 

where G na = n M G^ a represent specific components of 
the color field strength tensor (a being transverse). The 
zero momentum limit of this correlator is the gluonic pole 
matrix element mentioned above. It is the support of 
^g(x,x — xi) that we are after and specifically the model 
independent proof that it vanishes in the limit X\ — > 
for the case of the fragmentation correlators (|a;| > 1). 

The first step in all considerations is the observation 
that the k~ and fcj~ -integrations in the quark-quark and 
quark-quark-gluon correlators lead to light-front correla- 
tors, for which time-ordering is irrelevant. Therefore the 
matrix elements can be considered as forward matrix ele- 
ments of time-ordered products of operators. These rep- 
resent scattering amplitudes and their analytic structure 
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enables one to make statements about the support of the 
associated parton correlation functions. This can be done 
for quark-quark 27 j and multi-parton correlators [28j and 
for TMD as well as collinear correlators (29[. In this 
language the diagrams of Figs. [1] are just hadron-parton 
amplitudes, e.g. the quark-quark correlator Eq. (J3j) re- 
lated to the forward antiquark-hadron scattering ampli- 
tude A(k 2 ;s,u) (see Fig. [5]). Depending on the precise 
structure these are untruncated Green functions (time- 
ordered) or related to such Green functions via the LSZ 
formalism [30j |. 

The second step is the study of the analytic structure 
of an amplitude and in particular the singularities arising 
from cuts in the forward amplitudes. These are cuts in 
the (untruncated) legs, in particular the parton virtual- 
ly k 2 , and the Mandelstam invariants. The virtualities 
are conventionally placed just below the real axis, or the 
invariants are replaced by p 2 + ie, s + ie. The integra- 
tions over k~ and fcj" imply integrations over some of the 
invariants. At this point one must make the standard as- 
sumption that it is possible to use analyticity for QCD- 
amplitudes. We illustrate this step first for the standard 
quark-quark correlators. In that case one works with 
the four-point Green functions, shown in the middle of 
Fig. [21 where the cut amplitude depending on the particu- 
lar value of x gives the quark distributions (0 < x < 1) or 
the quark fragmentation functions (0 < z < 1 or i > 1) 
while analytic continuation to negative x describes the 
anti-quark distributions (— 1 < x < 0) and anti-quark 
fragmentation (-1 < 2 < or i < -1). 

The forward amplitude itself has singularities (cuts) 
for positive parton virtuality k 2 and in the Mandelstam 
variables, for which we choose s = (P — k) 2 and u = 
(P + k) 2 . For the forward amplitude, the invariants are 
constrained to s + u = 2 k 2 + 2 M 2 , where we will neglect 
the hadron masses (they don't play any essential role in 



s>o (u<0)j 



k < x < 1 



f u>o (s<0) 
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s > (u < 0) j 
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[ u>o (s<0) 



(b) antiquarks 



FIG. 3: The integration contours for the k~ integration with 
respect to the kinematic singularities in the (forward) anti- 
parton - hadron scattering amplitude for the case of (non- 
vanishing) distribution functions for quarks (a) and anti- 
quarks (b). 




(a) 



u - channel 



(b) 



FIG. 4: The additional invariants for the amplitude 
A(k 2 \ s, u\ 81, iti; k\, (k — ki) 2 ) relevant for gluonic pole matrix 
elements, (a) for the case s > and (b) for the case u > 0. 



our proof). These singularities then constitute the s-cut 
(for s > and u < 0) and the w-cut (for u > and 
s < 0). Using the expansion for k in Eq. @, one has 
k 2 = 2xk~ + k 2 or s = 2(x — l)k~ + k 2 and a similar 
expression for u. The transverse momenta, just as any of 
the parton or hadron masses, have little bearing on our 
results so we omit them in the expressions for k~ , 



s + ie 



2(x-l) 2(x + l) 



2.r 



(7) 



and one sees that for distribution functions the k~ inte- 
gration with respect to s, u and k 2 singularities follows 
the (dashed) contour in Fig. [3l The integration contours 
can be wrapped around the s and u-cuts for positive and 
negative .T-valucs respectively if \x\ < I, cuts that (in 
fcr) smoothly vanish when \x\ —> 1. Neither masses, nor 
transverse momenta matter and the support properties 
are valid for collinear and TMD PDFs. We get 

= 6(x) 6(1 - x) Disc [s] .4 + 6(-x) 6(1 + x) Disc [u] A. 

n (8) 

As discussed for instance in Ref. [31|, the case for frag- 
mentation is different since one in essence discusses the 
parton propagator for positive k 2 (sitting on the cut). 
For |x| > 1, one simply has 

A(a;) = 6(x - l)Disc [s] A + 6(-l - x)Disc [u] A. (9) 

In wrapping the integration around the s- or u-cut we 
have to assume convergence in the variable k~ (or k 2 ), 
or use subtracted relations. 

It is important to mention here that the integration of 
Eq. H]) and A(k 2 - s, u) over k T leads to ultraviolet diver- 
gences [2^, 28, 23| . However, this does not invalidate the 
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FIG. 5: The integration contours for the k^ integra- 
tion with respect to the singularities in the amplitude 
A(k 2 ; s, u\ si, m; k\ , (k — ki) 2 ) relevant for gluonic pole con- 
tributions. The figure shows for a given value of positive s 
(relevant for x > 0) how the fc-j" integration bypasses the cuts 
in si, ui and the parton virtualities in the limit xi — > +0. 
The cases < x < 1 (distributions) and x > 1 (fragmenta- 
tion) are shown in Figs (a) and (b) respectively. In the latter 
case only parton virtualities are relevant. 



assumption that the integral over k~ alone is sufficiently 
well behaved when k + and k T are fixed. Integration over 
k T as well as weighting with k T is anyway intimately 
linked with QCD evolution of TMDs 

We can extend this analyticity analysis to the multi- 
parton distribution and fragmentation functions in 
Eq. ([6]), by looking at the multi-parton amplitude 
A(k 2 ; s, u; Si,m; kf, (k — fci) 2 ) shown in Fig.|4j by study- 
ing the contours for the additional integrations. Defining 
the momentum of the additional parton as k\ as shown 
in Fig. [1] one retains the definitions and relations for s 
and u. For given positive s (s-channel, x > 0) or posi- 
tive u (u-channel, x < 0) one gets additional invariants 
Sl = (P =p k ± fci) 2 and u x = (P =F fci) 2 (cf- Figs H (a) 
and (b), respectively). Note that t\ = k 2 in both cases. 
Furthermore one has parton virtualities. Depending on 
if one is dealing with the s- or u-cut discontinuities, one 
has slightly different constraints for s± + ui, but for given 
values of s and u in the two cases of Fig. [4j one has cuts 
along si > (iti < 0) and u\ > (si < 0) as well as for 
positive parton virtualities. The relevant singularities for 
are found from 



■si + ie 



+ k~ = 



it 



2(n-(iTl)) 2(xiTl) 
k\ +ie (k- fci) 2 + ie 



2xi 



2{xi - x) 



(10) 



(with =p referring to s- and it-channel cuts, respectively). 
Again, parton or hadron masses as well as transverse 
momenta have little bearing as all they do is move the 
endpoints of the cuts. However, it is important to note 
that kiT in the complete expression for the numerator 



of Eq. (fTtj|) protects against the cut starting at zero in 
the zero mass limit. Depending on the value of xi, the 
integration contour in fcj~ bypasses the singularities en- 
countered in the complex plane in a particular way, which 
dictates the support properties of the quark-gluon-quark 
correlation functions. The denominators in Eq. (|10|l in 
the expressions relating fcf to si and u\ tell us that only 
whenxi G [re— 1,1] (for positive x) or x\ £ [— l,x+l] (for 
negative x) the singularities in s\ and ui are relevant. We 
study the case of the s-channel (x > 0). The it-channel 
is analogous. Looking at gluonic poles, we consider the 
limit X\ — » 0. For < x < 1, the value Xi = lies in 
the interval for which the Si and tti discontinuities can 
contribute. These are shown in Fig. [^a), now together 
with the singularities arising from the parton virtualities 
k\ and (k — ki) 2 . For the case x > 1 only these par- 
ton virtualities matter, shown in Fig.[5]^b). We then find 
(including for $g also the u-channcl contribution) in the 
limit x\ — > 



<S> G (x,x) 
A G (x,x) 



9(x)9{l - x) Discr ssi i.4 



+6(-x)6(l 
= 0, 



x)Bisc [u . Ul] A 1 (11) 
(12) 



where for Ag(x,x) the fcf integration can be wrapped 
around the k\ cut, which smoothly vanishes for xi — > +0. 
This is described by the arrow inside the branch cut in 
Figs. [5ja) and (b), indicating that it harmlessly recedes 
to infinity. Moreover it matches continuously to the case 
that xi < 0. Starting from x\ — > —0 one immediately 
would have obtained Eqs. (fTTj) and (fT2"|) , since the k\ cut 
is then along the negative fcf axis. This establishes the 
proof that gluonic pole matrix elements for fragmenta- 
tion correlators vanish. Similar results were obtained in 
the spectator model field theory in Ref. [lij], the spec- 
tator approach in Ref. fl7j | . and the general spectral ap- 
proach of Ref. [l8|. But now it has been demonstrated 
in a completely general way by assuming unitarity and 
analyticity properties to hold for QCD. Similar to the 
earlier discussion, we note that the integration of Eq. ((6|) 
and A{k 2 \ s, u\ si , tti ; k 2 , (k — fci ) 2 ) over ki T is ultraviolet 



divergent [2g, [28|, |29( . Again, the integral over fcf alone 
is well behaved when k^ and ki T are fixed with again full 
integration over k\ T corresponding to appropriate regu- 
larisation of A(k 2 ; s, u; Sj, iti; k 2 , (k — fci) 2 ) [33( and study 
of the QCD evolution [32j. 

In our last step, we show that our arguments for 
vanishing gluonic pole matrix elements hold for general 
multi-gluonic and even multi-partonic pole matrix ele- 
ments. Considering the analytic properties of general 
multi-gluonic pole matrix elements we can proceed in- 
ductively. For two gluons one simply extends the nesting 
of momenta k — k\ and fci by a nesting fc — fci — fc2, fci — k^ 
and fc2, which adds to the set (s, it, si, u\) two new invari- 
ants (s2, 112), without changing the behavior in the oth- 
ers. The gluonic pole matrix element /S.gg{x,x,x) thus 
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disappears as do all higher pole matrix elements. Since 
these higher pole matrix elements appear in the higher 



fc^-moments of the correlator A^(z,k T ) in Eq. (QJ, we 



conclude based on our very general assumptions of ana- 
lyticity for QCD amplitudes that this TMD correlator is 
universal and will be convoluted with the standard par- 
tonic cross sections. There is no proliferation of functions 
originating from the structure of the gauge links 14 [. 



This universality thus applies to all TMD fragmentation 
functions, T-even or T-odd, and for quark as well as for 
gluon PFFs. Most well-known are the T-odd ones for 
quarks such as the Collins function N. Il5i and the 
polarization fragmentation function D^ T [4j, |34j . These 
functions are simply allowed T-odd parts in the fragmen- 
tation correlator A(z,k T ) being a decay function. The 
corresponding T-odd TMD distribution functions, the 
Boer-Mulders function h^- and the Sivers function 
originate from the difference <£>I W +1 (x,p T ) — (x,p T ) 
of correlators with different gauge links and as a con- 
sequence will be convoluted with non-standard gluonic 
pole cross sections [l3|, 14 , [35{ . Our result also implies 



universality for the TMD fragmentation functions of glu- 
ons [361 ]; including for instance the T-even TMD frag- 



-L(S) 



mentation functions H 1 
ing distribution function 



link dependence |35l 137 



functions H^ 9 ^, however, is universal. In the case of 
these T-even functions, the non-trivial gauge link depen- 
dence only becomes visible in even fc T -moments involving 
contributions from T-even multi-gluonic pole matrix el- 
ements with an even number of gluons, all of which for 
fragmentation functions, however, will vanish. 
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, which just as the corrcspond- 
has a non-trivial gauge 
381 ] . The T-cven fragmentation 
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